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Using X-ray photoelectron spectroscopy, thermal desorption spectroscopy, and scanning tunneling 
microscopy we show that upon keV Xe^ irradiation of graphene on Ir(lll), Xe atoms are trapped 
under the graphene. Upon annealing, aggregation of Xe leads to graphene bulges and blisters. 

The efficient trapping is an unexpected and remarkable phenomenon, given the absence of chem¬ 
ical binding of Xe to Ir and to graphene, the weak interaction of a perfect graphene layer with 
Ir(lll), as well as the substantial damage to graphene due to irradiation. By combining molecular 
dynamics simulations and density functional theory calculations with our experiments, we uncover 
the mechanism of trapping. We describe ways to avoid blister formation during graphene growth, 
and also demonstrate how ion implantation can be used to intentionally create blisters without 
introducing damage to the graphene layer. Our approach may provide a pathway to synthesize new 
materials at a substrate - 2D material interface or to enable confined reactions at high pressures 
and temperatures. 

PACS numbers: 68.65.Pq, 61.80.Jh, 68.55.Ln, 68.55.J-, 68.37.Ef 


I. INTRODUCTION 

When graphene (Gr) adheres to a substrate, Gr blisters 
or bubbles can be formed by intentionally pressurizing a 
cavity in the substrat^®^, by agglomeration of atoms or 
molecules residing or being created in the space between 
Gr and the substrat^^^^, or even as a result of thermal 
processing. Such bubbles and blisters in Gr are of con¬ 
siderable interest in view of potential applications rang¬ 
ing from microscopic adaptive-focus lense^via fluid cells 
for transmission electron microscopy imagin^ISl or high 
pressure/ high temperature chemistr}® to surfaces with 
properties switchable by pressur^. Blister and bubble 
formation have also received considerable interest from a 
more fundamental point of view, since the strain associ¬ 
ated with them gives rise to huge pseudomagnetic fields 
in and thus to a quantization of the electronic 

structure into Landau levels. 

Bubble formation, blistering, and exfoliation in conse¬ 
quence of implantation of light elements like He, D, or H 
are phenomena that have been thoroughly investigated 
in the context of nuclear and fusion materialJ^. These 
phenomena were turned into useful technology by Bruel 
when applied to achieve microslicing of Si-wafers in sili¬ 
con on insulator technolog}®. Similarly, He+ irradiation 
of graphite in the basal plane orientation with a few ten 
keV ion energy has been found to result in gas bubble 
formation, blistering, and exfoliatior®!^. Using low en¬ 
ergy B and N ion implantation, substitutional doping of 
freestanding Gr was pioneered by Bangert et a/.® With 
similar low ion energies substitutional N implantation of 


Gr on SiG was accomplished by Telychko et al]^. For low 
energies and fiuences (negligible damage). Gun et alW^^ 
demonstrated that implanted Ar+ ions may come to rest 
under a 2D-layer strongly adhering to a metal surface and 
remain there upon annealing. In a recent commenfP^ to 
a paper by Herbig et al^, we pointed out that Xe+ irra¬ 
diation of Gr on Ir(lll) is accompanied by Xe trapping 
at the interface. 


Here we demonstrate that even a defective Gr layer 
on Ir(lll) effectively traps Xe supplied by ion irradia¬ 
tion up to very high temperatures. Upon annealing the 
trapped Xe aggregates cause flat bulges in Gr that trans¬ 
form eventually into highly pressurized blisters. Using 
a combined experimental and theoretical approach, we 
uncover that efficient Xe trapping is favored by colli- 
sional (reduced ion reflection due to Gr cover) and chem¬ 
ical (strong Gr edge bonds with Ir) effects. By proper 
noble gas pre-implantation and thermal processing, we 
grow noble gas filled Gr blisters with tunable size, areal 
density, and without introducing ion irradiation damage. 
The observed phenomena are expected to take place for a 
variety of substrates, ion species, and 2D-layers. We envi¬ 
sion that implantation of atoms with low solubility into a 
substrate covered with Gr (or another 2D-layer) followed 
by annealing could be used for stimulating confined re¬ 
actions in the space between Gr and the substrate, and 
as 2D-layers are transparent to light, for confined photo¬ 
chemical reactions. Such an approach could also be used 
for reactive growth of interfacial layers in between the 
substrate and a 2D-layer. 
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II. METHODS 

The experiments were conducted in ultra high vac¬ 
uum setups in Cologne [scanning tunneling microscopy 
(STM), thermal desorption spectroscopy (TDS)] and 
at the 1311 beamline of Max IV Laboratory in Lund 
[X-ray photoelectron spectroscopy (XPS)]. Cleaning of 
Ir(lll) was conducted by cycles of Xe+ (Cologne) or 
Ar+ (Lund), irradiation at room temperature (r.t.), O 2 
treatment at 1120 K and annealing to 1520 K (Cologne) 
or 1300K (Lund). A well-oriented and closed Gr layer 
was grown by one cycle of r.t. C 2 H 4 adsorption until 
saturation, subsequent thermal decomposition at 1300 K 
followed by Gr growth using 100 L ethylene at 1170 KP^. 
For the experiments represented by Fig. and Fig. S 6 
of the Supplemental Material (SMp^both temperatures 
were set to 1200 K. The ion flux during Xe+ irradiation 
at 300 K was controlled through adjusting the sample 
current that itself was calibrated for each ion energy by 
a Faraday cup in the Cologne setup. Ion incidence was 
~ 30° with respect to the surface normal for the XPS ex¬ 
periments (Lund) and at normal incidence for all other 
experiments (Cologne). Sputtering yields and ion ranges 
deviate by less than 10 % for these two angles of inci¬ 
dence. Annealing was conducted for 300 s at 1000 K and 
120 s at 1300 K, if not specified otherwise. The resulting 
sample morphology was imaged by STM at r.t. with a 
typical sample bias Ug ~ —IV and tunneling current 
It ~ InA. All XP-spectra were collected at r.t. in 
normal emission using a photon energy of 1000 eV for 
Xe3d and lr4s with a total energy resolution of light 
and analyzer better than 400 meV. After subtraction of 
a polynomial background the spectra were normalized to 
the height of the lr4s peak. TDS was conducted with a 
heating rate of 5K/s. 

The classical molecular dynamics (MD) simulations 
were run with the PARCAS cod^^ including only the 
atomic interactions, as for low ion energies electronic 
stopping can be neglected to good approximation. For 
each energy and system, 96 to 300 ions were shot, each 
separately, in perpendicular direction towards the system 
at a randomly chosen impact point. To model the sub¬ 
strate, we used Pt {Z = 78, mass 195 u) instead of an Ir 
substrate {Z = 77, mass 192 u) because a well-established 
interaction model exists for Pt-Pt and Pt-CP^, unlike for 
Ir-lr and Ir-G. Due to the similarity in atomic masses, 
structure, and chemistry of these two species, only minor 
differences are expected in the MD simulation results. 
The simulation cell size and time were adapted to allow 
sufficient relaxation of the system prior to analysis. The 
cell consisted of 18360 to 45900 Pt atoms to which 2584 
G atoms were added for Gr/Pt(lll). The simulation 
time ranged from 0.3 ns to 0.5 ns. Heat dissipation at 
the system edges was modeled with the Berendsen ther- 
mostafP^. The distance between the Gr sheet and the 
substrate was set to 3.31 A in agreement with the exper¬ 
imentally determined height of 3.3 Before the final 
relaxation the substrate lattice constant was stretched by 


1 % in order to coincide the two subsystems. The Pt-Pt 
and Pt-C interactions were modeled by the potentials of 
Albe et those of C-C by the bond-order potential 

of Brenner et al^ with a repulsive parfP^ for small atom 
separations, and the ion interactions by the universal re¬ 
pulsive potentiaP^. 

The first-principles density functional theory (DFT) 
calculations were performed using the plane -wave-basis- 
set Vienna ab initio simulation packagj^^^. The projec¬ 
tor augmented wave approaclP^ was used to describe the 
core electrons and a non-local van der Waals functionaP^ 
to describe exchange and correlation. A plane wave ki¬ 
netic energy cutoff of 400 eV was found to converge en¬ 
ergy differences between different configurations within 
0.1 eV. The same accuracy was achieved with regard to 
the number of k-points ( 3 x 3 x 1 ) in the two-dimensional 
Brillouin zone. All structures were relaxed until atomic 
forces were below 0.002 eV/A. The calculations were car¬ 
ried out for a 200 -atom 10 x 10 Gr supercell on top of 
a 9 X 9 three-atomic-layer-thick Ir (lll) slab containing 
243 atoms, as in our previous worlP^^, or for the same 
Ir slab without Gr. 20 A of vacuum was added in the 
transverse direction to separate the periodic images of 
the slab. The simulation supercell used for the calcu¬ 
lations of the migration path/energy of Xe from under 
graphene is shown in SM^. In practice, the Xe atom 
was placed in several positions along the path, and the 
atom was allowed to move only perpendicular to Ir(lll) 
during geometry optimization. All other atoms were al¬ 
lowed to move without any constraints. 


III. RESULTS AND DISCUSSION 
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FIG. 1. (color online) XP-spectra of the lr4s, Xe3d3/2, and 
Xe3d5/2 core levels for (a) bare Ir(lll) and (b) Gr covered 
Ir(lll) after exposure to 0.1 MLE of 0.5keV Xe^ at 300 K and 
annealed to stepwise increasing temperatures, as indicated. 
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In the first set of XPS experiments we exposed a bare 
Ir(lll) sample to 0.5 keV Xe+ at 300 K and conducted 
successive annealing to 1000 K and 1300 K. The ion fiu- 
ence was 0.1 MLE, where 1 MLE is 1.57 x lO^^ions/m^, 
i.e., numerically identical to the surface atomic den¬ 
sity of Ir(lll). The bottom XP spectrum in Eig. [^a) 
displays the Ir 4s core level peak together with the 
Xe3d3/2/Xe3d5/2 doublet after irradiation at 300 K. The 
Xe3d5/2 binding energy (BE) of 669.7 eV is consistent 
with values measur ed for implanted or adsorbed Xe at 
other metalJ^SE®^. Since Xe does not adsorb on Ir(lll) 
at 300 K (see our calculations below and ref. |43|), the 
Xe 3d signal must be attributed to Xe implanted into 
the Ir sample. Upon annealing to 1000 K the Xe 3d 
signal diminishes and has vanished after annealing to 
1300 K [middle and top spectra in Eig. [^a)]. Consis¬ 
tent with si milar studies for noble gas irradiation of other 
metalJ^^EH, we explain these changes as follows: due to 
thermal excitation Xe is released from its capture sites 
inside the crystal, diffuses to the surface, and finally des¬ 
orbs to the vacuum, such that only little (1000 K) or no 
(1300 K) Xe is left in the crystal. 

Conducting precisely the same irradiation experiment, 
but for Ir(lll) fully covered by Gr [hereafter named 
Gr/Ir(lll)], yields dramatically different XP spectra as 
shown in Fig. [^b). Already after irradiation at 300 K 
[bottom spectrum of Fig. [^b)] the integrated Xe3d5/2 
intensity is an order of magnitude higher (compare Ta¬ 
ble [I|. All the additional intensity (and thus nearly all 
intensity) must be due to Xe trapped in between the 
Gr cover and the Ir substrate. Annealing to 1000 K and 
1300 K hardly changes the integrated Xe3d peak inten¬ 
sity. The interpretation is straightforward: trapped Xe 
did not escape to any significant extent from under Gr 
to the vacuum, despite the fact that the Gr cover was 
heavily damaged by the Xe+ irradiation. 

Table compiles the integrated Xe3d5/2 intensities 
after ion exposure and annealing for bare Ir(lll) and 
Gr/Ir(lll). The intensities are given in % of the Xe3d5/2 
intensity of an adsorbed, saturated Xe layer on bare 
Ir(lll). Assuming a Xe saturation coverage of 0.33 ML 
on and neglecting attenuation effects, division 

of the tabulated intensities by a factor of three yields the 
amounts of trapped Xe in ML with respect to Ir(lll). 
Thus, after annealing to 1000 K, 4.1 %ML Xe is estimated 
to be trapped, i.e., more than 40% of the incident ions. 
This value is a lower bound, since the Ir 4s photoelectron 
attenuation through the adsorbed Xe layer for Xe/Ir(lll) 
and the attenuation of the Xe 3d signal from deeper lay¬ 
ers in Ir are neglected. Both factors tend to increase our 
estimate. Irrespective of these uncertainties, it is obvi¬ 
ous that a substantial fraction of the incident Xe remains 
underneath Gr up to the highest annealing temperature 
of 1300 K. 

While XPS tells about the Xe remaining at the sample 
after irradiation, a complementary view is provided by 
TDS, which yields information about the Xe leaving the 
sample after irradiation through thermal activation. The 


TABLE 1. Integrated Xe3d5/2 intensities /bare and Igt of 
bare and Gr covered Ir(lll), respectively, after exposure to 
0.1 MLE 0.5 keV or 3keV Xe^ at 300 K as well as annealing 
to 1000 K and 1300 K. Also, the intensity ratio /cr/Zbare is 
specified for the ion energies and temperatures used. Inten¬ 
sities are given as % of the integrated Xe 3d5 /2 intensity of a 
saturated Xe layer on Ir(lll) and are calibrated to the Ir 4s 
peak height measured simultaneously. 


0.5 keV 3keV 



/bare 

/Gr 

/Cr/ /bare 

/bare 

/Gr 

/Cr/ /bare 

300 K 

1.3 

11.5 

9 

2.9 

5.9 

2 

1000 K 

0.4 

12.2 

34 

1.9 

6.0 

3 

1300 K 

< 0.2 

11.6 

> 60 

1.5 

8.0 

6 


solid red line in Fig. is the thermal desorption trace of 
Xe after 0.5 keV irradiation of bare Ir(lll) with similar 
parameters as used for the XPS experiment. It exhibits a 
double peak structure with peak desorption temperatures 
of 745 K and 880 K. Beyond 1000 K the desorption rate 
gradually increases up to the end of the heating ramp at 
1330 K. The dashed red line in Fig. is the thermal des¬ 
orption trace of Xe after irradiation of Gr/Ir(lll). It is a 
fiat featureless curve with its slope increasing slightly for 
temperatures beyond 1000 K. Despite an order of magni¬ 
tude more Xe in the Gr/Ir(lll) sample, the integrated 
amount of Xe desorbing is lower by about a factor of two. 
Therefore, our TDS data elucidate that only a small frac¬ 
tion of the total amount of Xe in the Gr/Ir(lll) sample 
desorbs up to 1300 K. 



FIG. 2. TDS of Xe (131 amu) after exposure of bare Ir(lll) 
and Gr/Ir(lll) to 0.1 MLE Xe^ at 300K. Solid gray line: 
0.5keV, bare Ir(lll); dashed gray line: 0.5keV, Gr/Ir(lll); 
solid black line: 3keV, bare Ir(lll); dashed black line: 3keV, 
Gr/Ir(lll). Heating rate 5K/s. For clarity, the 3keV data is 
shifted vertically and multiplied by a factor of 0.5. 

The experimental data are substantiated by our MD 
simulations (compare Table , for which we used a Pt 
instead of an Ir substrate (see Methods section). For 
0.5keV Xe hitting bare Pt(lll) at normal incidence, 
when the energetic phase of the collision has ended (0.5 ns 
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after the impact), 98% of the impinging Xe has been re¬ 
flected (i.e., returned to vacuum) and only 2 % of the pri¬ 
mary ions have been implanted into the Pt crystal (into 
a depth of less than three atomic layers). The same sim¬ 
ulations with Xe hitting Gr/Pt(lll) yield a strikingly 
different picture. Then, only 2% of the projectiles are 
reflected and 5% are implanted into Pt, but 93% of the 
Xe ions are stuck in between Gr and the Pt(lll) surface 
[see Fig. SI of the SM for depth distribution^^ . Appar¬ 
ently, the primary Xe ions penetrate the Gr sheet easily; 
much more easily than a Pt layer. The reason for this is 
twofold: (i) due to the much smaller nuclear charge of C 
{Z = 12) compared to Pt {Z = 78) the scattering cross 
section of Xe with C is much smaller than that with Pt 
(and Ir); (ii) even in a head-on collision with a G atom the 
Xe atom loses less than a third of its energy and continues 
to move towards the metallic substrate. There, the Xe 
projectiles effectively transfer energy to the substrate Pt 
atoms (mass ratio 131 u to 195 u), such that in case of mo¬ 
mentum reversal they no longer possess enough energy to 
pass through the Gr sheet and are instead trapped. Note 
also that when impinging the Gr sheet from below, the 
Xe projectiles have not only much less energy and thus a 
larger scattering cross section with C, but under nearly 
all circumstances off-normal or grazing direction, such 
that the transparency of the Gr sheet for ion passage is 
strongly reduced. The results of the MD simulations are 
fully consistent with the XPS observation that an order 
of magnitude more Xe remains within the sample if there 
is a Gr cover present (compare Fig. and Table [^. Still, 
they do not explain why the trapped Xe does not escape 
to the vacuum by diffusion through irradiation induced 
vacancies and vacancy islands in the Gr sheet. Note also 
that the MD simulations are conducted for a perfect Gr 
sheet, while in the experiment ion induced damage of 
successive impacts accumulates. 



FIG. 3. (color online) (a) System energy as a function of 
the reaction coordinate (0 corresponds to the Xe atom at the 
Gr edge position) for two different edge structures and the 
path of a Xe atom schematically depicted in the inset, (b) 
Structure of zig-zag edge (left) and zoomed structure with 
the Xe atom in the transition state (right), (c) Structure of 
the reconstructed Klein edge (left) and zoomed structure with 
the Xe atom in the transition state (right). 


TABLE II. Fraction of reflected, trapped, or implanted pri¬ 
mary Xe ions given in % as obtained by MD simulations > 95 
events per case. Also specified is the average depth of implan¬ 
tation fim measured from the topmost Pt layer. 



0.5 keV 


3keV 


pt(iii) 

Gr/Pt(lll) 

pt(iii) 

Gr/Pt(lll) 

reflected 

98 

2 

57 

6 

trapped 

- 

93 

- 

41 

implanted 

2 

5 

43 

53 

fim (nm) 

0.4 

0.2 

1 

0.9 


To understand why trapped Xe remains under a de¬ 
fective Gr layer, we carried out DFT calculations, as 
described in the Methods section. We first evaluated 
the adsorption energy of Xe atoms on bare Ir(lll) and 
Gr/Ir(lll) to be —0.21eV and —O.lTeV, respectively. 
These numbers are consistent with previous experimen¬ 
tal results for Xe adsorption on Pt(lll) (—0.27 eV]P^and 
Xe adsorption on graphite (—0.25eV). The adsorption 
energies also agree up to a factor of about two with more 


accurate calculation schemes using many-body pertur¬ 
bation theory for Xe adsorbed to graphite, to Gr, or to 
Ni(lll), see e.g. l47f[49l In any case, the DFT numbers 
calculated here indicate the absence of chemical bonding 
and are in magnitude more than a factor of 30 smaller 
than the barriers for the penetration of Xe atoms through 
the graphene edge, as shown below. Next we inserted a 
single Xe atom between the Gr sheet and Ir(lll). After 
full geometry optimization, taking the energies of iso¬ 
lated Gr/Ir(lll) and the Xe atom as the reference, we 
found that it costs 2.9 eV to place a Xe atom between Gr 
and Ir (see Fig. S2 of the SM for geometr}Ell) . This en¬ 
ergy penalty results from the deformation caused by the 
Xe atom in the Gr sheet (elastic energy) and the local 
delamination of the Gr sheet from its substrate (delam¬ 
ination energy). Note that escape of a trapped Xe into 
the Ir bulk is impossible due to the very high intersti¬ 
tial formation energy exceeding 10 eV. The most likely 
escape path for Xe is thus to go via Gr edges onto the 
free Ir terrace, where the Xe atom would desorb rapidly 
under all temperatures of concern, as schematically de- 
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picted in the inset of Fig. [^a). Due to Gr sputtering, 
such edges are formed in large number during irradiation 
and annealing, but Gr dangling bonds are saturated due 
to the strong interaction of Gr with the metal giving rise 
to bending of the sh eet tow ards the metal surface already 
at room temperatur j^^ l ^^ i The energy of the system as a 
function of the coordinate of the Xe atom along the path 
is presented in Fig. |^a) for two particular edge struc¬ 
tures. First, the zig-zag edge that terminates vacancy 
islands in the Gr sheet [compare Fig. S5 (c) and caption]. 
Additionally, the reconstructed Klein edge is considered 
here, as such edge structures develop for freestanding Gr 
under sputtering condition^. For both structures an ac¬ 
tivation energy of more than 8eV is obtained. Fig. [^b) 
and (c) show the atomic structure of Gr edges binding 
to Ir(lll), and the close-ups present the Gr edge struc¬ 
ture and the Xe atom in the transition state. For passage 
through the zig-zag edge [Fig.j^b)] 2 C-Ir bonds must be 
broken temporarily, while for passing through the recon¬ 
structed Klein edge [Fig. 3(c)] two G-C bonds need to be 
broken in the C pentagons bound to the substrate. We 
also calculated the barrier for escape through an open 
corner in a Gr vacancy island edge and obtain a similar 
energy barrier (see Fig. S3 of the SIVP^. For an efficient 
release of Xe at 1300 K a much lower barrier, below 4 eV, 
would be necessary. Although computational cost pre¬ 
vents us from exploring the entire space of possible edge 
structures that may arise during irradiation and anneal¬ 
ing, it is clear from our calculations that the passage of a 
Xe atom from under Gr involves breaking of strong Gr-Ir 
or Gr-Gr bonds. Hence, very little Xe is released even 
through a highly defective Gr sheet with a large num¬ 
ber of smaller and larger vacancy islands as well as for 
temperatures up to 1300 K. As the escape of Xe from un¬ 
der Gr is impaired by Gr edge binding, the fraction of 
trapped Xe calculated in the MD simulations is indeed a 
good first estimate for the experimentally observed trap¬ 
ping efficiency. 

The real space view of the noble gas trapping under the 
Gr sheet is presented in the STM topographs of Fig. 
After room temperature irradiation the sample displays 
a rather dense distribution of small scale protrusions of 
0.3 nm - 0.5 nm height [compare Fig. |^a)]. Experiments 
with lower fluence show unambiguously that the protru¬ 
sions are deformations in the Gr sheet, which itself is 
defective but globally intacfP^. Based on our XPS and 
TDS experiments, as well as the MD simulations and 
DFT computations discussed above, we must interpret 
these protrusions to result to a large extent from Xe 
atoms, which give rise to Gr’s local deformation. Ad- 
ditionall y, I r adatoms pushed onto the Ir(lll) surface by 
collisionP^ contribute to the formation of protrusions in 
the Gr sheet. This statement is backed-up by the inset 
of Fig. [^a), where Ir adatom clusters are visible after 
irradiation of bare Ir(lll) using the same ion beam pa¬ 
rameters. 

After annealing the sample to 1000 K the morphology 
is dramatically changed [compare Fig. [^b)]. The num¬ 


ber of protrusions is strongly reduced and two types of 
elevations are visible: first, flat bulges with a well defined 
height around 0.35 nm, and second, larger, smoothly 
curved blisters with a maximum height up to 0.8 nm 
[height profile in Fig. [^b)]. Also visible are dark ar¬ 
eas of few nm size and point-like dark spots with typical 
separations of 2-3 nm, as highlighted by the inset. These 
features are vacancy islands within the Gr sheet, either 
extended nm-sized areas or small point-like spots pinned 
to preferential binding sites within the moire formed of 
Gr with the Ir(lll) substrate (compare reference l38|) . 
The vacancy islands result from aggregation of vacancies 
created in Gr due to sputtering. 

Upon annealing to 1300 K the morphology coarsened 
and is now dominated by few, large blisters with lat¬ 
eral dimensions of the order of 10 nm and heights around 
1.5 nm, as shown in Fig. [^c). Also some shallower blis¬ 
ters are present. The point-like vacancy clusters have 
nearly vanished. Large vacancy islands are present in 
quite low density. One of those is visible in the upper 
right corner of Fig. |^c), between two large blisters. Ap¬ 
parently the small vacancy clusters detached during an¬ 
nealing from their preferential binding sites and aggre¬ 
gated to large vacancy islands, consistent with the find¬ 
ings of Standop et al^. It is striking that the vacancy 
island visible in Fig.|^c) is bound by large blisters. This 
observation underlines the strong binding of Gr to the 
Ir(lll) substrate and the impermeability of the Gr edge 
for Xe, as found in our DFT calculations. A rather regu¬ 
lar moire corrugation with a 2.5 nm pitch indicates that 
the Gr lattice has restored itself completely [compare in¬ 
set of Fig. [^c)]. 

We tentatively interpret the evolution from atomic 
scale protrusions to large blisters as follows. After 
room temperature irradiation, a mixture of Xe atoms, 
Ir adatoms, and detached G atoms is present under the 
defective Gr cover. Upon annealing to 1000 K the Ir 
adatoms and adatom clusters recombine with Ir surface 
vacancies restoring a flat surfac^^, while the detached G 
atoms are largely reincorporated into the Gr sheet. The 
agglomeration of the remaining Xe atoms (and of the 
implanted atoms diffusing from the Ir bulk to the sur¬ 
face) gives rise to monolayer Xe platelets that cause flat 
bulges. Their formation is driven by minimization of de¬ 
formation and delamination energy within Gr: disperse 
single Xe atoms cause a much larger deformation of Gr 
and much more delamination of Gr from the substrate 
than one flat aggregate containing all Xe atoms. The as¬ 
sumption of a monolayer thick Xe aggregate is consistent 
with the uniform apparent STM bulge height of about 
0.35 nm, similar to the Xe van der Waals diameter of 
0.43 nm. Since the Gr sheet presses the Xe together, we 
tentatively assume that also the lateral spacing of the Xe 
atoms is close to that of solid Xe. This could be real¬ 
ized, e.g., W a (v^S X \/3)-Xe structure with a coverage 
of 0.33 MlP^. When the flat bulges grow larger, a shape 
transformation to smoothly curved blisters takes place. 
While the delamination energy grows with the square 
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FIG. 4. (color online) STM topographs of Gr/Ir(lll) and associated height profiles after exposure to 0.1 MLE of 0.5keV Xe^ 
(a) at 300 K, (b) annealed to 1000 K, and (c) annealed to 1300 K. The height prohles below the topographs are taken along the 
lines indicated. Insets: (a) bare Ir(lll) after exposure to ion irradiation with the same parameters, (b), (c) are zooms of the 
corresponding topograph. Image sizes are 90 nm x 90 nm, inset sizes are (a) 24 nm x 24 nm and (b), (c) 15 nm x 15 nm. 


of the linear dimension of the flat bulge, its strain and 
bending energy - primarily localized along the rim of the 
bulge - grows only proportional to it. Therefore, at some 
point a transition to a more three-dimensional shape of 
the Xe cavity takes place, lowering the delamination en¬ 
ergy for a given amount of Xe. After annealing to 1300K 
the measured blister height and radius (approximating 
the blister as spherically symmetric) with the calculated 
adhesion energy of Gr to Ir(lll) (50meV per C atonP^ 
allows one to estimat^^ pressures of the order of IGPa 
for the Xe in the blister. It should be noted that a similar 
coarsening of gas cavities during annealing, as described 
here, is also observed for bubbles and blisters after noble 
gas implantation into bulk material^^. 

To explore how trapping depends on ion energy, we 
conducted the same set of experiments with an ion en¬ 
ergy of 3keV. Despite the increased Gr sputtering and 
damaging (the yield for C atom detachment from the 
Gr sheet increases from 1.0 to 1.8 according to our MD 
simulations) the general picture is unchanged: the Gr 
sheet strongly reduces ion reflection and enables efficient 
Xe trapping upon irradiation at 300 K. To analyze differ¬ 
ences to the 0.5 keV situation, we first consider changes 
for the irradiation of bare Ir(lll). Compared to 0.5 keV, 
for 3keV Xe+ into bare metal substantially more Xe is 
implanted (more than 40 % of the impinging Xe atoms) 
into a larger average depth of 1 nm, as apparent from 
the MD results shown in Table nn and consistent with 
the corresponding XPS data provided as Fig. S4 of the 
SlVpIl (compare also ref. [24]). Consequently, substantially 


more Xe is released from bare Ir during heating as vis¬ 
ible in the corresponding Xe desorption trace in Fig. 
However, even after heating to 1300 K some Xe is left in 
the Ir crystal, presumably due to incorporation of Xe in 
stable bulk vacancy clusters formed in consequence of the 
more violent Xe impacts. With the Gr cover the implan¬ 
tation into Ir does not change much, but trapping under 
the Gr sheet becomes frequent (> 40% of the impinging 
Xe), thereby strongly suppressing reflection. As for the 
0.5 keV case, the presence of the Gr cover gives rise to 
a substantial enhancement of the Xe 3d XPS intensity 
after 300 K irradiation, which must be attributed to Xe 
trapping (compare Fig. S4 of the SIVP^ and Table |^. 
Xe release is largely prevented by the Gr cover, as obvi¬ 
ous from the comparison of bare Ir(lll) and Gr/Ir(lll) 
Xe desorption traces after 3keV ion irradiation shown 
in Fig. 1^ In consequence, by annealing, the amount of 
trapped Xe continuously increases by Xe streaming to the 
interface from its implantation sites in the Ir bulk. The 
STM images display a markedly larger Gr removal result¬ 
ing from the higher sputtering yield of 3 keV Xe ions, but 
otherwise a similar morphological evolution from atomic 
scale protrusions to large blisters (compare Fig. S5 of 
the SIVPI). 

Our experimental analysis, as well as the simulations 
and calculations presented here indicate that the trapped 
Xe fully accounts for the bulges and blisters observed. 
Though we cannot rule out the presence of additional in¬ 
terfacial C, the association of bulge formation with inter¬ 
facial C (Ref. [25]) appears not to be plausible. More likely. 
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upon annealing the interfacial C reincorporates into the 
Gr sheet. The trapping yields presented in reference!^ in 
Figs. 5 and 6 should therefore be interpreted as Xe trap¬ 
ping efficiencies^. The observation of bulges and blisters 
after Ar+ and Ne+ irradiation of Gr/Ir(lll) in reference 
[25] indicates efficient trapping under the Gr sheet also for 
these noble gases. Only for He+ exposure with 0.3 keV 
flat bulges and blisters are almost absent, even after an 
ion fluence of 0.3 MLE, indicating that trapping does not 
take place. Presumably, due to the small van der Waals 
diameter (0.28nm for He compared to 0.43nm for Xe), 
He escapes from under Gr without the need to break Gr 
edge bonds. Indeed, our DFT calculations of He diffu¬ 
sion from under the Gr sheet with the same path as for 
Xe showed that the barrier for He escape is considerably 
smaller, on the order of 3-4 eV. Sputtering with He+, 
preferably even at elevated temperature, may therefore 
be a method to avoid noble gas trapping during ion irra¬ 
diation of Gr on a substrate. 

In our view, the phenomenon of trapping after irra¬ 
diation of 2D-materials is of general relevance. It does 
not depend sensitively on ion beam parameters (energy, 
angle of incidence, species, charge state), as long as the 
ions are able to penetrate the 2D-layer. It should also 
take place for other combinations of substrate and 2D- 
layer materials, as long as (i) the implanted species does 
not react with substrate and 2D-layer, (ii) the solubility 
of the implanted species in the substrate is low, and (hi) 
the 2D-layer interacts strongly with its defects and edges 
to the substrat^^. 

Our investigations indicate that high temperature 2D- 
layer growth on an improperly degassed substrate with 
a higher than equilibrium gas concentration may result 
in inferior quality of the 2D-layer through the formation 
of gas cavities (blisters, bubbles, wrinkles) at the inter¬ 
face. For example, Gr grown on a metal substrate that 
was sputter-cleaned with keV noble gas ions and not an¬ 
nealed sufficiently long well above the Gr growth temper¬ 
ature, is likely to display gas filled blisters after growth. 
An example for blister formation during Gr growth on a 
sputter-cleaned and insufficiently annealed Ir(lll) sub¬ 
strate is shown as Fig. S6 of the 

In the experiments described up to now, blister forma¬ 
tion is a by-product of ion irradiation and always involves 
damaging of the Gr sheet. However, it is also possible to 
separate ion irradiation and blister formation, enabling 
us to create blisters under a perfect Gr sheet with num¬ 
ber density and size controllable by thermal processing. 
To reach this goal we first implanted Xe at 300 K into 
bare Ir(lll), using sufficiently high and varied ion en¬ 
ergies to ensure high gas release temperatures and gas 
incorporation with a broad depth distribution. To re¬ 
move surface damage, the sample was flash annealed to 
1200 K and subsequently a Gr sheet was grown at 1200 K. 
Fig.j^a) displays the resulting morphology characterized 
by smooth circular bulges of subatomic height due to gas 
bubbles in IiH]^ 

as well as irregularly shaped (often trian¬ 
gular) blisters that formed due to gas release from the Ir 


crystal during Gr growth at 1200 K. The Xe still present 
in the subsurface bubbles can be released and trapped 
underneath the perfect Gr sheet by additional annealing 
to 1500 K, as shown in Fig. [^b). The bubbles within 
the Ir crystal vanished, while the volume of the blisters 
increased by a factor of five. 

We envision that similar procedures of pre¬ 
implantation and annealing could be used to create 
pressurized nanoreactors with one reagent supplied 
through release of implanted species from the substrate. 
Also the growth of interfacial layers, e.g., insulators, 
could be accomplished by providing a reagent from the 
bulk after 2D-layer growth or transfer and optional 
intercalation. 


IV. CONCLUSIONS 

In conclusion, we have determined the mechanisms of 
atom trapping underneath Gr on Ir(lll) during Xe+ ion 
irradiation. It relies on the combination of two effects. 
First, for an ion energy up to a few keV a Gr layer 
on Ir(lll) effectively decreases ion reflection in favor of 
ion trapping in between Gr and its substrate. The ef¬ 
ficient trapping results from the easy penetration of Xe 
through Gr (small scattering cross section) and efficient 
energy loss to the metallic substrate, making escape after 
momentum-reversal unlikely. Second, the irradiation in¬ 
duced Gr edges bind strongly to the Ir(lll) substrate 
making escape of Xe through out-diffusion at the Gr 
edge an unlikely process, even at annealing temperatures 
of 1300 K. Driven by elastic and delamination energy of 
Gr, upon annealing the trapped Xe first forms mono- 
layer thick Xe islands under fiat bulges in Gr and, when 
growing larger at higher temperatures, three-dimensional 
pressurized blisters in the Gr sheet. The bulge and blister 
formation upon ion irradiation must be assumed to take 
place for a broad range of 2D-materials which adhere 
to the substrate. We demonstrated that ion beam im¬ 
plantation methods in combination with proper thermal 
processing may be a pathway to enable reactions in pres¬ 
surized nanocavities or to grow homogeneous interfacial 
layers between a perfect 2D material and its substrate. 
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FIG. 5. (color online) (a) STM topograph of Gr/Ir(lll) after successive exposure to 0.2 MLE of 5keV, 4keV, 3keV, and 2keV 
Xe^ (accumulated exposure 0.8 MLE) at 300 K, subsequent annealing to 1200 K, and eventually growth of a complete layer of 
Gr at 1200 K. (b) sample after additional flash annealing to 1500 K. Image sizes are 250 nm x 250 nm. 
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